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i.  INTROmiCTION 


The  generation  of  sound  by  a  ^ght  source  was  suggested  more  than  one 
hundred  years  ago  by  Alexander  Sraham  Bell.^  It  has  received  more 
attention  recently  because  of  new  developments  in  high  power  lasers,  which 
make  this  acoustic  generation  mechanism  more  suitable  for  practical 
applications. 

Westervelt  and  Larson  have  shown^  ^at  ver.y  directive  sound  beams  of 
low  frequency  can  be  achieved  with  a  thermoacoustic  array  (or  optoacoustic 
antenna).  A  thermoacoustic  array  is  an  acoustic  line  source  obtained  by 
shining  a  laser  into  a  medium.  Optical  absorption  in  the  medium  induces  an 
exponential  tapering  of  the  line  source  along  the  axis  of  penetration  of 
the  laser  beam.  The  principle  of  operation  of  these  arrays  is  that 
amplitude  modulation  of  the  intensity  of  the  laser  beam,  directed  toward 
the  absorptive  medium  (usually  water),  induces  a  periodic  heating  of  the 
medium  and  therefore  a  fluctuation  of  density.  This  in  turn  generates  an 
acoustic  wave  in  the  medium.  The  acoustic  frequency  is  thus  equal  1  o  the 
modulation  frequency  of  the  laser  intensity. 

The  first  experimental  verifications  of  thermoacoustic  radiation  were 
4  5  6 

made  by  Muir  et  al .  ’  and  Culbertson  and,  as  expected,  the  efficiency  of 
the  conversion  of  electromagnetic  energy  into  acoustic  energy  was  extremely 
low.  However,  recent  studies  by  Bunkin  et  al.  ’  have  shown  that  the 
amplitude  of  the  acoustic  signal  can  be  increased  by  moving  the  laser  beam 
at  high  velocities  through  the  water.  High  speed  motion  of  the  beam  can 
easily  be  achieved  using  a  rotating  mirror;  therefore  a  mechanical 
transducer  moving  in  the  medium  is  not  required. 

This  report  describes  a  simple  time  domain  method  which  can  be  used  to 
compute  the  pressure  waveform  radiated  by  a  moving  thermoacoustic  array 
(MTA).  In  Section  II,  we  shall  derive  an  expression  for  the  acoustic 
response  observed  in  a  medium  when  the  optoacoustic  source  is  an  ideal 
impulse.  In  this  case,  the  acoustic  response  will  be  called  the  impulse 


response  of  a  thermoacoustic  array  (TA).  It  is  assumed  that  the  light 
source  is  in  the  form  of  a  narrow  beam  striking  the  medium  at  normal 
incidence.  In  Section  III,  the  impulse  response  of  a  TA  will  be  used  to 
compute  the  pressure  waveform  radiated  by  a  source  moving  at  any  velocity 
through  an  optically  absorbing  medium. 

There  are  three  main  reasons  for  choosing  the  impulse  response 
approach  rather  than  a  frequency  response  approach:  (1)  it  gives 

analytical  results  valid  in  both  the  nearfield  and  farfield  of  the  source, 

(2)  the  theory  does  not  break  down  for  transonic  motion  of  the  source,  and 

(3)  time  information  is  presented  without  need  for  a  transform. 

In  Section  IV  the  time  domain  method  is  compared  to  a  previous  theory 

g 

developed  by  Lyamshev  and  Sedov  for  the  special  case  of  farfield  radiation 
with  analytical  predictions  based  on  the  impulse  response  approach.  It 
will  be  shown  that  the  two  theories  are  essentially  in  agreement. 

The  main  conclusions  of  this  study  will  be  surmarized  in  the  last 
section. 


II.  IMPULSE  RESPONSE  OF  A  THERMOACOUSTIC  ARRAY 


In  this  section  an  analytical  expression  is  derived  for  the  impulse 
response  of  a  thermoacoustic  array.  The  first  step  is  to  solve  the  wave 
equation  and  derive  an  expression  in  integral  form  for  the  acoustic 
pressure  in  terms  of  the  optical  intensity  of  the  source.  The  impulse 
response  is  then  obtained  by  letting  the  source  function  be  replaced  by  a 
delta  function. 


A.  Analytical  Approach 


The  wave  equation  for  a  viscous  medium  containing  heat  sources  is 
derived  in  Appendix  A.  For  an  inviscid,  non  heat  conducting  fluid,  the 
linearized  wave  equation  takes  the  form 


JL  _ia 

cp  3t 


in 


where 

p  is  the  acoustic  pressure, 
c  is  the  small  signal  sound  speed, 

q=-  V.l(x,y,z,t)  is  energy  per  unit  volume  and  per  unit  time  added  to 
the  medium, 

I  is  the  laser  intensity, 

8  is  the  logarithmic  coefficient  of  thermal  expansion  of  the  medium, 
and 

cp  is  the  specific  heat  of  the  medium  at  constant  pressure. 

The  solution  to  Eq.  (1)  for  a  pressure  release  air-water  interface  is 
obtained  by  integrating  the  Green's  function  of  the  problem  over  the  source 
region. 
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Here  r  and  r'  are  distances  from  the  source  to  the  receiver  and  are  shown 
in  Fig.  1. 

Equation  (2)  gives  the  acoustic  pressure  p  observed  at  the  receiver  at 

any  point  in  the  field.  The  first  integral  in  Eq.  (2)  represents  the 

contribution  of  the  thermoacoustic  source  located  in  the  water,  and 

therefore  the  volume  of  integration  is  the  volume  V  .  of  the  column  of 

water 

water  being  irradiated  by  the  laser  beam.  The  second  term  in  Eq.  (2) 
represents  the  effect  of  the  image  sources  located  above  the  interface,  and 
therefore  the  volume  of  integration  is  Vair,  the  mirror  image  of  Vwater 
(see  Fig.  1.).  The  minus  sign  is  due  to  the  compliant  nature  of  the 
interface.  These  two  contributions  are  not  identical  because  the  time  lags 
r/c  and  r'/c  are  different  for  the  volumes  of  integration.  For  a  uniform 
intensity  distribution  across  the  beam,  the  energy  q  reduces  to 

q  =  -  -§7  [A  I(t)  e'aZ]  ,  (3) 

-CtZ 

where  A  is  the  optical  transmissivity  from  air  to  water,  and  e  is  the 
exponential  tapering  with  depth  associated  with  the  absorption  of  light  in 
water. 

Assuming  that  the  beam  diameter  is  always  very  small  compared  to  the 
characteristic  wavelength  X  of  the  acoustic  radiation,  we  rewrite  Eq.  (2) 
in  the  following  form: 
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where  S=fTa  is  the  beam  cross  section.  For  a  laser  pulse  whose  intensity 
is  modulated  at  frequency  f,  the  characteristic  acoustic  wavelength  is 
X  =  c/f ;  for  an  unmodulated  laser  pulse,  X=CTp,  where  xp  is  the  pulse 
duration. 


Combining  Eqs.  (3)  and  (4),  we  obtain  a  general  expression  for  the 
acoustic  pressure  in  terms  of  the  laser  intensity: 


P 


A  8  S  a 

_Tlrcp 


A  6S  a 

ATTCp 


(5) 


Equation  (1)  shows  that  the  optoacoustic  source  strength  (forcing 
function)  is  proportional  to  the  time  derivative  of  the  laser  intensity. 
It  is  therefore  convenient  to  define  the  impulse  response  h(t)  by 

P(t)  =  h(t)  *-|i  ,  (6) 

where  the  asterisk  denotes  convolution.  The  function  h(t)  is  therefore  the 
acoustic  response  to  an  impulse  dl/dt.  Consequently,  h(t)  can  be  regarded 
as  either  the  optoacoustic  impulse  response  or  the  light  step  response. 
From  Eq.  (5)  we  find 


h(t) 


A8  Sa 
4  mCp 


dz  -  mirror  image 


(7) 


Shifting  the  integration  variable  from  z  to  r,  we  rewrite  Eq.  (7)  as 


^  A  8  S  etc 
h(t)  = 

P 


mirror  image 


(8) 
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Equation  (8)  provides  a  means  for  evaluating  the  impulse  response  of  a 
thermoacoustic  array  once  the  relationship  between  the  depth  z  and  radius 
of  observation  r  or  r‘  is  known.  Assuming  that  the  laser  beam  strikes  the 
water  at  normal  incidence,  the  array  may  be  divided  into  four  regions. 
These  are  indicated  in  Fig.  2.  For  each  region  the  relationship  between  z 
and  r  or  r'  is  given  below. 


-region  1  (mirror  image): 


<  z  <0 


-  <  r 1  <  oo 

o  - 


?  2  2 

z  =  -rQ  cos0Q  +Vr  -  rQ  sin  0q 


dz 
dr 1 

-region  2: 


r  1 


2  2  2 
r  -  r  sin£0„ 
o  o 


0  <  z  <  r„  cos0„  ;  r^  Sin9^  i  r  <  rr 


o  o 


o  o 


2  2  2 

z  =  -r  cos0„  +  Vr  -  r  sin  0 
oo  oo 


dz 

37 


-r 


r  -  s1"  90 


-region  3: 


ro  cos0Q  <  z  <  2ro  cos0Q  ;  rQ  sin0o  <  r  <  rQ 


z  =  rQ  cos0Q  +  Vr2  -  r2  sin20o 


dz 

dr 


2  2  2 

r  -  rL  sin^9 
o  o 


•region  4: 


z  >2r„  cos0 


r  <  r  <  oo 


time  delay  i-At  taken  by  the  source  to  travel  from  the  origin  (t=0)  to 
the  ith  position.  This  gives  the  total  acoustic  pressure  Py(t)  radiated 
by  the  thermoacoustic  array  during  its  motion: 

N  MC 

PT(t)  =  £  Pi (t-i -At)A t  . 

1=0 

Combining  Eqs.  (18)  and  (19)  finally  gives 

N  ,  or 

p  (t )  =  £M'(i*At)  h.(t-i»At)At  . 

i=0 

Equation  (20)  shows  that  the  total  pressure  received  at  the  observation 
point  is  a  convolution  type  summation  in  the  time  domain,  between  the 
"unsteady"  (changing  shape  with  time)  impulse  response  h^ (t)  and  the 
optoacoustic  source  strength  dl/dt.  Note  that  in  the  case  of  a 
stationary  array,  the  shape  of  impulse  response  remains  constant  and  the 
numerical  program  performs  a  true  convolution  between  the  impulse 
response  and  the  optoacoustic  source  strength. 

This  time  domain  approach  is  perfectly  suitable  for  a  numerical 
computation,  and  a  program  to  perform  the  computation  has  been  written  on 
the  ARL:UT  CYBER  computer. 

B.  Farfield  Criterion 


The  farfield  criterion  has  been  a  source  of  confusion  in  the 
literature  on  thermoacoustic  arrays.  It  seems  necessary  to  state  clearly 
the  difference  between  the  geometric  farfield  (Fresnel  approximation), 
the  acoustic  farfield,  and  a  farfield  criterion  based  upon  the  source 
dimensions  (Fraunhofer  approximation),  for  both  stationary  and  moving 
sources,  with  or  without  intensity  modulation  of  the  laser  pulse.  The 
geometric  farfield  depends  upon  the  dimension(s)  of  the  source  relative 
to  the  distance  of  observation.  The  acoustic  farfield  depends  upon  the 
distance  of  observation  relative  to  a  typical  acoustic  wavelength.  The 
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The  second  step  is  to  compute  for  each  of  these  positions  the 

elementary  acoustic  pressure  p^(t)  radiated  by  the  array  when  it  was  at 
the  position  shown  on  the  ith  picture,  that  is  to  say,  after  a  time  delay 
i*At.  The  laser  source  intensity  at  that  instant  is  simply  I  ( i  *  A  t )  and 
the  impulse  response  h.(t)  can  be  evaluated  from  Eqs.  ( 10 ) - ( 14 ) . 
However,  it  is  very  important  to  realize  that  the  impulse  response  is  a 
function  of  the  position  of  the  array  because  the  distance  rQ  and  the 

angle  90  between  the  source  and  the  receiver  are  time  dependent  in  the 
case  of  a  moving  source.  Let  us  denote  by  r.  and  9,  the  values  of  r  and 

9q  shown  on  the  ith  picture  corresponding  to  time  i *At  after  the 

beginning  of  the  motion  of  the  source.  It  can  very  easily  be  shown  that 


r . 

l 


-  2ro(v.i*At)  cos§0  +  (v.i*At)^ 


(16) 


and 


(17) 


A.  A 

where  r  and  9  refer  to  values  at  the  origin  of  the  motion  of  the  source 
(t=0).  The  impulse  response  h-(t)  is  therefore  given  by  the  set  of 
Eqs.  ( 10 ) - ( 14 )  where  the  coordinates  rQ  and  9Q  have  been  replaced  by  r. 
and  9i  ,  using  Eqs.  (16)  and  (17). 


The  elementary  acoustic  pressure  p.(t)  is  then  computed  as  follows. 
Equation  (6)  shows  that  the  optoacoustic  source  strength  is  proportional 
to  the  time  derivative  of  the  laser  intensity.  Therefore  during  the  ith 


picture,  the  optoacoustic  strength  is 


t=i *At 


or, 


in  a  shorthand 


notation,  1 1 ( i * At ) .  The  ith  picture  therefore  shows  an  array  excited  by 

an  impulse  of  strength  I ' ( i -A  t ) ,  and  then  p.(t)  is  given  by 


Pi  (t)  =  I '  ( i-A  t)h.  (t) 


(18) 


The  last  step  is  then  to  add  all  the  elementary  acoustic  pressures 
p.(t)  rad.ated  during  the  motion  of  the  source,  taking  into  account  the 
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III.  MOVING  THERMOACOUSTIC  ARRAY 


It  has  been  shown  in  the  previous  section  that  the  acoustic  impulse 
response  of  a  thermoacoustic  array  is  described  by  Eqs.  (10)-(14).  In  this 
section  we  seek  a  solution  for  the  pressure  field  radiated  by  a  moving 
thermoacoustic  array  (MTA)  in  terms  of  its  impulse  response.  It  is  assumed 
in  what  follows  that  the  observation  point  is  on  the  plane  of  motion  of  the 
source. 

A.  Time  Domain  Analysis 

The  pressure  radiated  by  an  MTA  can  be  constructed  in  a  three-step 
procedure.  First,  the  motion  of  the  source  can  be  decomposed  in  time.  As 
in  a  movie,  a  "picture"  can  be  taken  at  constant  time  intervals  At, 
showing  the  laser  beam  at  different  positions  along  its  path  during  the 
laser  pulse  duration  Tp.  Second,  for  each  of  these  positions  (defined  by 
the  subscript  i)  the  impulse  response  h^ (t)  and  the  corresponding 
elementary  pressure  response  p^(t)  may  be  computed  numerically  from 
Eqs.  (10)-(14).  Third,  the  total  pressure  p^(t)  received  at  the  hydrophone 
can  then  be  obtained  by  adding  all  of  the  elementary  acoustic  responses 
p^t)  with  the  suitable  time  delays  corresponding  to  the  motion  of  the 
source.  Note  that  this  approach  is  also  suitable  for  a  source  moving  with 
a  nonuniform  velocity  along  its  path;  however,  in  the  following  analysis, 
it  will  be  assumed  for  simplicity  that  the  source  is  moving  at  a  constant 
velocity  v. 

The  first  step  in  the  analysis  is  a  discretization  in  time  and  space 
of  the  motion  of  the  source  so  that  the  problem  can  be  analyzed 
numerically  on  a  computer.  Given  that  xp  is  the  laser  pulse  duration 
and  At  the  time  increment,  the  motion  of  the  source  is  represented  by  N 
pictures  taken  every  At  seconds,  such  that  N=xp/At  .  For  a  source 
moving  recti  1  inear ly  at  constant  velocity  v,  i •  A t  represents  the  time 
taken  for  the  array  to  move  from  the  origin  to  the  os  it  ion  shown  on  the 
ith  picture. 
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kw  w : v  v*j  ?? 


dh 

HT 


t=t  . 
min 


=  0 
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-  « 
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with  h(t)  defined  by  Eq.  (11),  leads  to  tanh(ay  m.n)  =  l/ay  ^n,  i.e.. 


4 

J 


ay  .  =1.2 

min 


where 


‘min  V^ctmin^  "  ro  sin2®o 


Figure  4  shows  the  evolution  of  the  impulse  response  function  h(t)  in 
the  r -t  plane. 
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strongest  disturbance  is  therefore  the  one  which  suffers  the  least 
spherical  spreading  loss.  Consequently  the  strongest  disturbance  is 

originated  by  the  elementary  source  on  the  array  which  is  the  closest  to 
the  receiver.  Since  the  minimum  distance  between  the  source  and  the 
receiver  is  rQsin90,  the  impulse  response  h(t)  exhibits  a  peak  at  time 
rQsin90/c  and  then  decreases  gradually  until  t=r  /c. 

It  is  also  interesting  to  note  that  the  first  effect  (exponential 
decay  in  the  source  strength  along  the  z  axis)  dominates  the  second  effect 
(spherical  spreading  loss)  when  the  observation  point  is  located  in  the 
farfield.  Reciprocally,  the  second  effect  becomes  predominant  when  the 
observation  point  is  situated  in  the  nearfield.  This  means  that  in  the 
farfield  most  of  the  acoustic  response  occurs  at  t=rQ/c,  while  in  the 
nearfield  most  of  it  is  expected  to  occur  at  t=rQsin9o/c.  To  show  this, 
let  us  recall  that  the  observation  point  is  said  to  be  in  the  geometrical 
farfield  when  the  vertical  diffraction  loss  of  the  array  can  be  neglected. 
This  will  be  discussed  in  more  detail  in  Section  III.B.  This  farfield 
approximation  is  in  fact  equivalent  to  neglecting  the  difference  in 
spherical  spreading  between  all  the  elementary  sources  of  the  array.  The 
remaining  mechanism  determining  the  shape  of  the  impulse  response  h(t)  is 
thus,  in  the  geometric  farfield,  the  exponential  tapering  along  the  z  axis. 

Reciprocally,  in  the  nearfield,  the  difference  of  path  from  the 
elementary  sources  of  the  array  to  the  observation  point  can  be  quite 
different  from  one  elementary  source  to  another.  The  spherical  spreading 
loss  associated  with  this  difference  of  path  is,  in  the  nearfield,  always 
more  important  than  the  exponential  shading  along  the  z  axis,  and  therefore 
dictates  the  shape  of  the  impulse  response  h(t). 

This  approach  may  be  used  to  define  a  geometrical  nearfield  criterion. 
When  both  spherical  spreading  and  exponential  shading  are  of  the  same 

importance,  the  impulse  response  exhibits  two  peaks  of  approximately  equal 
magnitude  at  time  t=rQsin9o/c  and  t0=r0/c  (see  Fig.  3  for  the  case  T  =  3). 
Thus  h(t)  passes  through  a  minimum  value  at  a  time  t  ^n,  such  that 

tl<tmin<to  evaluatin9 
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and  I  and  K  are  the  modified  Bessel  functions  of  first  and  second  kind  of 
n  n 

order  n. 


B.  Discussion 

The  following  section  is  a  detailed  analysis  of  the  behavior  of  the 
impulse  response  h(t)  when  r0sin90/c<t<r0/c,  whose  aim  is  to  provide  some 
understanding  of  the  nearfield  of  a  thermoacoustic  array.  The  shape  of 
h(t)  can  change  drastically  from  a  monotonically  decreasing  function  with  a 
maximum  at  t=rQsin90/c,  to  a  monotonically  increasing  function  having  a 
maximum  at  t=rQ/c.  This  curious  behavior  can  be  explained  by  looking  at 
the  two  different  effects  determining  the  shape  of  the  impulse 
response:  the  exponential  shading  e"  az  along  the  axis  of  the  array,  which 
is  due  to  the  absorption  of  light  in  the  medium,  and  the  spherical 
spreading  associated  with  each  wavelet  radiated  by  the  thermoacoustic 
array. 


For  simplicity  let  us  first  assume  that  the  exponential  decay  along 
the  depth  of  the  TA  is  much  more  important  than  the  spherical  spreading. 
Let  z^,  zB,  Zq,  and  zQ  be  the  depths  of  the  elementary  sources  located  at 
points  A,  B,  C,  and  D  on  the  array  as  shown  in  Fig.  2.  It  can  be  shown 
that,  for  any  positive  a  , 


azr 


-  az. 


-  az, 


azr 


+  e 


<  e 


+  e 


Therefore  the  strength  of  the  acoustic  disturbance  coming  from  points  B  and 
C  is  always  less  important  than  the  strength  of  the  one  originating  at 
points  A  and  D.  Since  the  latter  arrives  at  the  receiving  point  at  time 
rQ/c ,  the  impulse  response  h(t)  is  expected  to  increase  gradually  from  time 
r0sin90/c  to  time  rQ/c. 

If  we  assume  now  that  the  spherical  spreading  loss  is  much  more 
important  than  the  exponential  decay  of  the  source  strength  along  the 
array,  all  elementary  sources  have  approximative^  the  same  strength.  The 
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FIGURE  3 

IMPULSE  RESPONSE  OF  A  THERMOACOUSTIC  ARRAY 


with 


„  _  A  8S  ac 
K  -TTE^- 


lo  -r  • 

r^sinO 
.  oo 


r  =  arQ  cosQ0  ,  and 

_  r  /2  12 

y  “  c  vt  ^  1  " 

Figure  3  shows  how  the  shape  of  the  impulse  response  h(t)  changes 
as  the  nondimensional  parameter  T  =arocos90  is  varied  from  0.02  to  100. 
(  r  is  the  ratio  of  the  depth  of  the  receiver  to  the  effective  length  a-1 
of  the  thermoacoustic  array.)  It  is  assumed  that  the  laser  is  shined  into 
fresh  water  with  c  =  1486  m/sec. 

The  impulse  response  h(t)  may  also  be  expressed  in  terms  of 
modified  Bessel  functions.  The  relationship  between  hyperbolic  functions 
and  modified  Bessel  functions  is  given  in  Ref.  10,  p.  443.  It  can  easily 
be  shown  that 

f  0  for  t  <  t. 


K  e"r  (rosin9o)_1 


2 it  Ka  e"r  P'1/2  I  ,  /9(p) 


for  t  =  t. 


for  t,  <  t  <  t 
1  o 


K  e"2F  (ro)_1 


for  t  =  t. 


_l  k  a  s  i  nh  {  T  )  p'1/2  K1/2(p)  for  *  >  ^ 


where 


p  =  ay  =  avfct)2  -  sin20o 
13 


(13) 


.  -  ^  WV~*  V"»V*  V'.  V  W  \  x  V  *_  •*.  -. 


h(t)  =  K 


-2ar  cosQ 
0  0 


5.  t  >  rQ/c 

The  disturbances  emitted  by  regions  1  and  4  of  the  TA  will  arrive 
by  pairs  at  the  observation  point.  However,  in  this  case  the  boundary 
condition  at  the  air-water  interface  (pressure  release)  implies  that 
signals  coming  from  region  1  (mirror  image)  are  inverted.  In  this  time 
region  the  only  nonzero  contribution  to  the  impulse  response  is  given  by 
the  last  two  integrals  in  Eq.  (9).  These  may  be  evaluated  to  yield 


h(t)  =  .  ■  .  ■>  ==r  exp 

V^ct)2  -  r2  sin29Q  l 


b(ro 


cosOq  +V  (ct)‘ 


-  ro  slAo). 


-  exp  a  (-rQ  cos0Q  +  ^/(ct)2  -  r2  sin2  -  9Q))_  J 
which  can  be  rewritten  in  simpler  form  by  using  the  hyperbolic  sine, 

-2K  exp  [-evict)2  -  rj;  sinZ90]  . 


h(t)  = 


^ct ) 2  -  r2  sin20o 


s i nh (  arQ  cos0Q)  .  (14) 


The  impulse  response  of  a  thermoacoustic  array  as  a  function  of 
time  is  therefore  given  by  the  set  of  Eqs.  (10)-(14). 


0 

for 

t 

Eq. 
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K  e" 
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r1 
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t 
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Eq. 
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-r  -l 

y 

cosh(ay  ) 
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tl 

A 
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A 

c+ 
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Eq. 

(12) 

K  e~ 

2F  (r  )-X 
v  o' 
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t 

=  »o 

Eq. 

(13) 

-2K 

-1  -ay 
y  e 

sinh(  r  ) 

for 

t 

>  tn 

Eq. 

(14) 

3. 


<  t  < 


r  sinQ 
o _ o 

c 


r 

o 
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By  referring  to  Fig.  2,  it  is  obvious  that  the  first  wavelet  that  will 
reach  the  receiver  emanates  from  point  E.  A  short  instant  later,  wavelets 
coming  from  points  B  and  C  will  arrive  simultaneously  at  the  receiver.  The 
progression  continues  with  wavelets  eventually  arriving  from  points 
A  and  D.  These  are  the  wavelets  received  during  the  interval 
rQsin90/c<t<r0/c. 


The  disturbances  coming  from  regions  2  and  3  of  the  TA,  then,  will  add 
coherently  by  pairs  at  the  observation  point.  The  only  nonzero 
contribution  to  the  impulse  response  is  given  by  the  first  two  integrals  in 
Eq.  (9),  which  may  be  evaluated  to  give 

h(t)  =  /  o  K  o  ?  jexp  [-a (r  cosQ  -  ct )2  -  r2  sin20Q  ) 

y(ct)2  -  r2  sin20o  l  L  V  J 

+  exp  £-a'(r0  cos0Q  +^/(ct)2  -  r2  sinZ0Q  ) 

This  expression  can  be  simplified  by  introducing  the  hyperbolic  cosine. 


-  a  r  cos0  _ , 

h(t)  =  ~7"~2  2  2  cosh[a  V{ct)Z  ~  ro  sin20o]  '  (12) 

V^ct)  -  r0  sin  9o 


4.  t  =  rQ/c 

There  are  two  points  on  the  array  which  are  located  at  a  distance 
rQ  from  the  receiver.  These  two  points  are  denoted  A  and  D  on  Fig.  2.  The 
impulse  response  at  t=rQ/c  is  thus  the  sum  of  both  wavelets  coming  from 
points  A  and  D.  The  pressure  release  boundary  condition  at  the  free 
surface  of  the  liquid  implies  that  the  acoustic  contribution  from  point  A 
is  zero.  Therefore  the  impulse  response  at  t=r  /c  is  simply  given  by  the 
wavelet  coming  from  D,  that  is  to  say  from  the  point  of  the  array  located 
at  a  depth  z=2rocos0Q.  Equation  (7)  yields  an  expression  for  the  impulse 
at  t=r  /c. 


t  <  ro  Sin0o/c, 
t  =  rQ  sin90/c, 
rQ  sin0Q/c  <  t  <  rQ/ c, 
t  =  rQ/c,  and 
t  >  rQ/c. 

Using  Eq.  (9)  it  is  possible  to  find  the  impulse  response  h(t)  for  the 
five  time  regimes. 

1.  t  <  rQ  sin0o/c 

The  impulse  response  reduces  to 

h(t)  =  0  ,  (10) 

which  states  that  the  acoustic  emission  from  the  array  has  not  yet  had  time 
to  arrive  at  the  source.  The  first  disturbance  emitted  by  the  TA  has  to 
travel  a  distance  rQ$in0o  before  reaching  the  observation  point.  Therefore 
the  impulse  response  h(t)  is  expected  to  be  identically  zero  for  any  time 

t<rosin9o/c' 

2.  t  =  rQ  sin0Q/c 

Since  ^sinOg  represents  the  shortest  distance  between  the  source  and  the 
receiver  (see  Fig.  2),  the  particular  instant  t=rQsin0o/c  represents  the 
time  of  arrival  of  the  first  wavelet  at  the  receiver.  This  wavelet  was 
originated  at  a  distance  rQs i n0Q  from  the  receiver,  that  is  to  say,  at  a 
depth  z=rQcos90  on  the  array.  Consequently,  the  impulse  response  is, 
simply,  from  Eq.  (7) , 

-“Vos0o 

h(t)  =  '  r0  sin80-  •  HD 
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z  =  r0  cosd0  +vr  -  rQ  sin  9Q 


dz  _  _ r _ 

W  "/“p  2  T2 

-  ro  sin  9o 


Expanding  Eq.  (8)  in  separate  integrals  over  the  four  regions  leads  to 


h(t)  =  K 


r  sinQ„ 
o  o 


['  “(r° 


cos9oV^^\)  /2-fdr2 

Vr  -  rQ  Sin  9 


r~9.tI  exp  |-a  ^rQ  cosQQ\/r^ 


-r0  5 'A,) 


rosin9o 


/r2  -  r2  sin290 


'  “(ro 


ro 


cos0oVr  -  r 


o  si  A)] 


'  -  r0  s1  A 


/r2  .  r2  .  2  ,(9 
r  r0  sin^Q 

0 


where 


A  SSctc 


The  four  integrals  on  Eq.  (9)  represent  the  contributions  to  the 
acoustic  pressure  from  the  four  array  regions  defined  previously.  However, 
it  is  more  physically  meaningful  to  distinguish  the  array  response  by 
regions  in  time  rather  than  in  space.  The  limits  of  integration  indicate 
that  the  analysis  of  the  impulse  response  should  be  divided  into  five  time 
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FIGURE  2 

GEOMETRY  OF  THE  ARRAY 


source  dimensions  approximation  depends  upon  the  dimensions  of  the  source 
relative  to  a  typical  acoustic  wavelength. 

1.  Geometric  Farfield 


a.  Stationary  Source 


It  can  be  seen  from  Fig.  5  that,  for  an  array  of  effective 
length  L=cT'*‘,  the  radial  distances  to  the  receiver  from  the  array  depths 
z=0  and  z=L  are  related  according  to 

2-J/2 


rL  =  ro 


1  -  2  ft)  C0S9°  +  ftj 

.11 


The  Fresnel  approximation  for  the  geometric  farfield  is  given  by 


(21) 


r<< 1 

o 


(22) 


If  this  approximation  is  met,  Eq.  (21)  may  be  rewritten  as 


rL  =  rQ  -  L  cos90 


(23) 


We  may  define  a  time  delay  as 


r  -  r,  L  cos  9 
o  L  o 


(24) 


This  time  delay  is  characteristic  of  the  vertical  diffraction  across  the 
effective  length  of*  of  the  array  projected  onto  the  direction  of 
observation.  It  can  be  seen  from  Eq.  (22)  that  large  values  of  the 
nondimensional  parameter  r=r0cos90/L  indicate  a  geometric  farfield 
behavior.  It  is  therefore  a  measure  of  the  geometric  nearfield  of  a 
stationary  TA. 
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b.  Moving  Source 


If  the  source  is  moving  at  velocity  v  in  the  +x  direction 


during  a  time  duration  t  it  travels  a  distance  £=vTp. 


Now  the  source 


may  be  thought  of  as  rectangular  with  dimensions  Lx&  .  In  order  to  be  in 
the  geometric  farfield,  both  dimensions  must  be  small  relative  to  r  . 
Hence  l  will  play  a  role  similar  to  that  of  L  in  the  previous  discussion. 
The  Fresnel  approximation  for  the  geometric  farfield  is  then 


or,  if  M=v/c, 


Similarly  we  can  define  a  time  delay  characteristic  of  the 
horizontal  diffraction  due  to  the  motion  of  the  source  (see  Fig.  5), 


ro  “  r* 


where  r?  is  the  distance  to  the  receiver  from  the  array  point  at  z=0  and 
x=£.  For  £«rQ,  -^reduces  to 


£sin9 


=  M  t  sinQ 
P  o 


Equations  (22)  and  (26)  define  the  geometric  farfield  criterion.  These 
equations  are  valid  for  both  modulated  and  unmodulated  source 
intensities. 

2.  Acoustic  Farfield 

The  acoustic  farfield  criterion  is  a  measure  of  the  ratio  of 
the  distance  of  observation  to  a  typical  wavelength  emitted  by  the 
source.  The  acoustic  farfield  approximation  is  therefore  completely 
independent  of  the  geometric  farfield  approximation. 


a.  Stationary  Source 

When  the  distance  of  observation  rQ  is  much  bigger  than 
the  wavelength  X  emitted  by  the  source,  the  observer  is  said  to  be  in  the 
(acoustic)  farfield.  This  occurs  when 


kr0  >  >  1 


where  k  is  the  wave  number  of  the  lowest  frequency  component  present  in 
the  signal  (most  restrictive  case). 


imposed  by  Eq.  (29)  becomes 


For  an  unmodulated  pulse,  k=2TT/cTp,  so  the  condition 


t  <  <—  (3 

pc  ' 

For  a  monochromatic  modulation  of  a  signal  at  a  frequency 
fQ  (period  TQ) ,  k=2  Tr/cTQ,  so  that  the  farfield  condition  can  be 

expressed  as 


b.  Moving  Source 


T  <  < 

o  c 


The  previous  section  has  just  shown  that  the  acoustic 

farfield  is  determined  by  the  relative  importance  of  the  typical  time 

scale  of  the  pulse  to  the  retarded  time  of  observation  r  /c.  However, 

o 

when  a  source  is  in  motion,  it  is  well  known  that  the  time  scale  will 
change  in  accordance  with  the  Doppler  shift.  The  acoustic  farfield 
criterion  can  therefore  be  rewritten  as 


k^r_  >  >  l 


(3 


where  is  the  Doppler  shifted  wave  number  of  the  lowest  frequency 

component  in  the  signal.  In  the  case  of  a  point  source  moving  at 
velocity  v=Mc  in  the  +x  direction,  is  related  to  k  by 


k 


d 


k 

1 1-M  sin0Q  ' 


(33) 


Using  Eq.  (33)  as  an  approximation  for  the  wave  number  of  the 
thermoacoustic  array,  the  acoustic  farfield  criterion  for  an  unmodulated 
pulse  of  duration  xp  becomes 


>  >  C  T 


1-M  sinQ, 


(34) 


For  a  pulse  modulated  every  TQ  seconds,  the  acoustic  farfield  occurs  when 


T 


o 


1-M  sin9Q 


(35) 


Note  that  Eqs.  (34)  and  (35)  are  always  satisfied  in  the  Cferenkov  direction 
(sin9Q=M"'*')  so  that  for  a  source  moving  at  velocity  c/sin9Q,  the  acoustic 
nearfield  is  virtually  nonexistent.  However,  it  must  be  borne  in  mind  that 
Eqs.  (34)  and  (35)  are  only  crude  approximations  to  the  real  physical 
problem--the  Doppler  shift  introduced  in  these  expressions  is  valid  only 
for  a  point  source.  Moreover,  the  equations  assume  that  the  angle  of 
observation  remains  constant  during  the  motion  of  the  source.  This 
assumption  is  obviously  violated  at  large  Mach  numbers.  The  exact  Doppler 
shift  of  a  thermoacoustic  array  moving  at  any  velocity  (including 
transonic)  has  been  investigated  in  more  detail  and  will  be  presented  in 
a  separate  report. 


3.  Source  Dimensions 


In  Section  IV,  the  results  obtained  with  the  impulse  response 

technique  described  in  Sections  II  and  III  will  be  compared  with  some 

g 

results  obtained  in  a  previous  study.  This  comparison  will  be  easier  to 
make  if  we  introduce  the  Fraunhofer  approximation.  It  is  important  to 
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realize  that  this  approximation  is  independent  of  the  farfield 
approximations  described  previously.  The  Fraunhofer  approximation 
depends  only  upon  the  dimensions  of  the  source  relative  to  the  acoustic 
wavelength  radiated  by  the  source. 

a.  Stationary  Source 

Assume  that  the  source  is  small  in  terms  of  an  acoustic 
wavelength,  i.e.,  that 


kL  <  <  1  , 


(36) 


where  L  is  the  length  of  the  source  and  k  is  the  wave  number  of  the 

hiqhest  frequency  component  in  the  signal  (most  restrictive  case).  This 

11 

is  the  Fraunhofer  approximation. 

.  For  an  unmodulated  pulse,  k=2*/c'rp,  so  that  by  combining  Eq.  (24)  and 
Eq.  (36)  we  obtain 


V*  Tl 


(37) 


•  For  an  intensity  modulated  signal  of  period  TQ,  Eq.  (36)  can  be 
rewritten  as 


T  >>  t  .  (38) 

o  L 

Equations  (37)  and  (38)  represent  the  small  (stationary)  source 
approx imat i on  for  both  unmodulated  and  modulated  cases. 

For  a  source  that  is  large  in  terms  of  a  wavelength  (long 
narrow  cylinder). 


kl  >  >  1 


(39) 


i 
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where  k  is  the  wave  number  of  the  lowest  frequency  component  of  the 
signal.  For  both  modulated  and  unmodulated  cases,  k=2Tr/CTp,  and 
Eq.  (39)  can  be  rewritten  as 


<  < 


tL 


(40) 


b.  Moving  Source 

Here  again  the  time  scale  change  (Doppler  shift)  will 
affect  the  characteristic  wave  number  in  Eq.  (34)  and  Eq.  (37)  (refer  to 
Section  III . B . 2 . b) .  As  in  the  case  of  the  geometric  farfield 
approximation,  the  motion  of  the  source  introduces  a  second  dimension  to 
the  source.  Therefore  the  source  is  small  in  terms  of  a  wavelength  if 
both  the  following  equations  are  satisfied. 

kdL<<  1  (41) 


and 


kd*«  1  .  (42) 

The  small  source  approximation  for  an  unmodulated  pulse  of  duration 
becomes 


P.d 


>>  T 


L’ 


and,  for  a  modulated  pulse  of  period  T  ,  becomes 


>  >  T 


L’ 


(43) 


(44) 


where  the  subscript  d  indicates  a  Doppler  shifted  quantity.  Similarly, 
if  the  source  is  long  in  terms  of  a  wavelength. 


and 


Ji 


<  <  T, 


(46) 


Equation  (45)  indicates  the  characteristic  of  a  long  source  in  the  z 
direction,  while  Eq.  (46)  indicates  the  property  of  a  long  source  in  the 
x  direction.  Both  Eqs.  (45)  and  (46)  apply  for  either  modulated  or 
unmodulated  laser  pulse. 


C.  Results 


In  this  section  we  present  the  results  obtained  using  the  impulse 
response  technique  to  predict  the  pressure  waveform  radiated  by  a 
thermoacoustic  array.  We  first  discuss  the  laser  pulse  shape  which  has 
been  chosen  for  the  numerical  predictions.  We  then  restrict  our 
attention  to  the  results  obtained  when  the  source  is  stationary  and  the 
intensity  unmodulated.  Finally  we  show  the  theoretical  predictions 
obtained  in  some  realistic  configurations  of  a  moving  and  modulated  laser 
pulse. 


1.  Laser  Pulse 


5 

Muir  et  al .  have  shown  that  a  modulated  laser  intensity  may  be 
conveniently  described  by 

I (t )  =  IQ(t)  sin2^)  =  IQ(t)  [l/2  -  1/2  coswt]  ,  (47) 

where  IQ(t)  is  the  laser  pulse  envelope,  and  oj  is  the  angular  modulation 
frequency,  that  is,  the  acoustic  frequency  for  a  stationary  source. 

The  laser  pulse  envelope  IQ(t)  was  determined  from  the  actual 
laser  used  in  the  experimental  study  of  thermoacoustic  arrays  at  ARL : UT 
facilities.  It  was  found  that  IQ(t)  is  well  represented  by 

yt)  =  A  (y  exp 
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with  A=10.8  and  B=5.0.  The  peak  amplitude  occurs  at  t=ip/B,  where  is 
the  laser  pulse  duration. 

2.  Stationary  Unmodulated  Laser  Pulse 

The  impulse  response  approach  described  previously  was  used  to 
predict  the  pressure  waveform  radiated  by  a  stationary  thermoacoustic 
array.  The  modulation  frequency  was  set  to  zero  in  order  to  simplify  the 
interpretation  of  the  results.  However,  it  should  be  understood  that  the 
computer  program  is  not  restricted  to  unmodulated  pulses. 

Figure  6  shows  the  evolution  of  the  predicted  pressure  waveform 

received  in  fresh  water  at  the  observation  point  as  the  nondimensional 

parameter  r  changes  from  200  to  0.2,  i.e.,  from  geometric  farfield  to 

nearfield  values.  The  position  of  the  receiving  point  is  kept  constant 

(r  =  4  m  and  9  =60°)  so  that  the  variation  of  T  corresponds  in  fact  to 
u  u  1 

the  variation  of  the  effective  length  (L=  a"A)  of  the  array,  from  1  cm  to 
10  m  (short  array  to  long  array).  The  pulse  duration  was  set  to  250  usee 
so  that  the  observation  point  was  always  located  in  the  acoustic  farfield 
(kro=67). 

In  the  farfield  case  shown  in  Fig.  6(a),  the  pressure  waveform  is 

found  to  occur  over  a  narrow  time  period  and  have  a  sharp  peak  at  the 
onset,  followed  by  a  negative  signal.  The  sharp  peak  in  the  acoustic 
response  observed  in  Fig.  6(a)  is  due  to  the  fact  that  the  laser  pulse 

described  by  Eq.  (48)  has  a  discontinuous  time  derivative  at  t=0.  The 
inverted  shape  of  the  laser  pulse  which  appears  for  all  values  of  T  in 
Figs.  6(a)-6(f)  after  the  retarded  time  rQ/c  can  be  explained  by  the 
effect  of  the  pressure  release  interface  between  air  and  water.  As  will 
be  shown  in  Section  IV. A  of  this  report,  this  farfield  waveform  is 
directly  proportional  to  the  second  time  derivative  of  the  laser  pulse 

intensity. 


As  shown  in  Figs.  6(b)-6(f),  the  pressure  waveform  broadens  and 
changes  substantial ly  in  shape  as  the  situation  approaches  the  geometric 


r  =  200 


r  =  4o 


FIGURE  6 

PRESSURE  WAVEFORM  AS  A  FUNCTION  OF  T  =  arQcos  6>0 
FOR  A  STATIONARY  UNMODULATED  LASER  PULSE 

(rQ  =  4  m;  0Q  =  60°;  rp  =  250  H sec) 
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nearfield. 


In  the  extreme  nearfield  represented  by  Fig.  6(f),  the 

pressure  waveform  is  seen  to  begin  at  time  t=rosin9o/c  rather  than  the 
time  t=rQ/c,  at  which  the  pressure  signal  begins  in  the  farfield.  This 
result  agrees  with  the  analysis  presented  in  Section  II. B.  It  is  also 
interesting  to  note  that  in  the  extreme  geometric  nearfield  of  a  very 
long  and  narrow  TA  (Fig.  6(f)),  the  acoustic  response  occurring  at  the 
retarded  time  r  /c  will  be  exactly  an  inverted  laser  pulse.  This  effect 
is  also  discussed  in  Section  IV. B. 

3.  Moving  Modulated  Laser  Pulse 

Some  typical  pressure  waveforms  radiated  by  a  moving 
thermoacoustic  array  are  presented  in  Fig.  7.  These  waveforms  were 

obtained  by  computer  simulations,  using  the  impulse  response  approach 
described  previously.  These  predictions  were  made  for  realistic  values 
of  the  important  parameters:  rQ  =  4  m,  90=75°,  a  =  15  m”\  fQ  =  7  kHz, 
Tp  =  1  msec,  and  T=15.5.  The  Mach  number  was  varied  from  0  to  1.5. 

Table  I,  which  shows  the  farfield  measures  as  a  function  of  the 

Mach  number,  indicates  that  the  observation  point  was  always  in  the 

acoustic  farfield  (kcjro>>^'  However,  it  was  in  the  geometric  nearfield 
for  M=1 . 04  and  M=1.5.  Table  I  also  shows  that  the  source  used  in 
generating  Fig.  7  was  always  long  in  both  the  x  and  z  directions  in  terms 
of  a  wavelength,  except  when  the  source  was  stationary. 

Figure  7(a)  shows  the  acoustic  responses  of  a  modulated,  stationary 
thermoacoustic  array.  There  is  a  periodic  nature  to  the  response  with  a 

roughly  exponential  decay  in  the  amplitude  corresponding  to  the  exponential 
decay  in  the  laser  intensity  (see  Eg.  (48)). 

Figures  7(b),  7(c),  and  7(d)  show  the  acoustic  response  of  a  moving 

modulated  thermoacoustic  array.  Note  that  these  figures  show  waveforms  of 

different  shapes  than  the  stationary  array.  Some  of  the  changes  caused  by 

the  motion  of  the  array  are  the  result  of  moving  into  the  geometric 
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FIGURE  7 

ACOUSTIC  RESPONSE  OF  A  MOVING  THERMOACOUSTIC  ARRAY 


TABLE  I 


FARFIELD  LIMIT  FOR  THE  THERMOACOUSTIC  WAVES  PLOTTED  IN  FIGURE  7 


nearfield.  Other  changes  are  caused  by  the  time  compression  which  occurs 
for  nonstationary  sources. 

One  advantage  of  the  time  domain  approach,  used  in  the  numerical 
computation  of  the  pressure  waveform,  is  that  the  Doppler  shift  is 

implicitly  taken  into  account  by  the  time  delays  of  the  pseudo-convolution. 

Therefore  the  Doppler  shift  appears  naturally  and  is  valid  not  only  for  a 
point  source,  but  also  for  an  extended  array  of  length  This  is 

discussed  thoroughly  by  Berthelot  and  Busch-Vishni ac,  who  show  that  the 
effect  of  having  an  extended  source  instead  of  a  point  source  is  to  reduce 

v 

the  pressure  amplitude  measured  in  the  Cerenkov  direction.  This  is  why  the 
computed  pressure  waveform  is  always  bounded  even  for  a  Mach  wave  (see 
Fig.  7(c)). 

The  case  of  a  thermoacoustic  array  moving  at  transonic  velocity 

deserves  special  attention  because  it  yields  the  maximum  peak  power  of  the 
optoacoustic  source.  A  detailed  analysis  of  the  Doppler  shift  associated 
with  a  transonic  source  may  be  found  in  Ref.  13.  The  difficulty  arising 
with  transonic  source  motion  is  the  assumption  that  the  laser  beam  diameter 
is  small  compared  to  an  acoustic  wavelength.  This  assumption  is  clearly 
violated  at  transonic  velocities  because  in  that  case  the  Doppler  shifted 
acoustic  wavelength  goes  to  zero.  It  is  therefore  essential  to  include  the 
effects  of  the  finite  width  of  the  laser  beam  on  the  pressure  waveform 
radiated  by  a  transonic  thermoacoustic  array.  A  report  describing  these 
finite  effects  is  being  prepared. 


In  order  to  analyze  Figs.  7(b)  and  7(d),  we  may,  however,  simplify  the 

problem  and  use  as  a  first  approximation  the  classical  Doppler  factor 
12 

formulation  without  any  finite  beamwidth  effects.  In  the  formulation, 
the  Doppler  factor  D  is  given  by 


D 


1  -  7  sin  9 
c  o 


(49) 


Equation  (49)  indicates  that  a  modulation  frequency  f  of  7  kHz  should 
give  a  Doppler  shifted  frequency  f d=f q/D  of  about  13.5  kHz  at  M=0.5  and 
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It  is  also  convenient  to  relate  the  energy  q,  added  to  the  medium,  to  the 

-f 

intensity  of  the  laser  beam  illuminating  the  medium.  Let  us  denote  by  <I> 
the  intensity  of  the  laser  beam  averaged  over  a  period  of  the  (very  rapid) 
optical  oscillation  of  the  laser  light.  < I >  is  called  the  Poynting  vector 
of  the  process.  Then  the  flux  of  intensity  V*<I>  represents  just  the  heat 
absorbed  by  the  medium  (per  unit  time  and  unit  volume).  Therefore  it 
follows  from  the  definition  of  q  that 


q  =  -  V  •  <  I  > 


(A. 17) 


In  practice  the  laser  beam  is  shining  in,  say,  the  +z  direction,  so  that 
9  I 

q  =  "  37  • 


Equations  (A. 16)  and  (A. 17)  are  the  starting  point  of  the  theory  presented 
in  this  report. 


It  is  interesting  to  note  that  Eq.  (A. 16)  has  exactly  the  same  form  as 
the  Wes tervelt  inhomogeneous  wave  equation ^  describing  parametric  end-fire 
radiation.  In  the  wave  equation  for  the  parametric  array,  however,  the 
source  term  represents  virtual  sources  generated  by  nonlinear  interaction 
in  the  medium,  whereas  in  Eq.  (A. 16)  the  source  term  represents  the  time 
and  spatial  derivative  of  the  intensity  distribution  in  the  laser  beam. 
There  is  nevertheless  a  close  mathematical  analogy  between  parametric 
arrays  and  thermoacoustic  arrays,  and  this  analogy  has  been  used  by 
Novikov,  Rudenko,  and  Timoshenko^  to  study  the  pressure  waveform  radiated 
by  a  transonic  thermoacoustic  array. 


(A. 10) 


-  ■  TT-'V  -1 


P  =  C 


'P  -  ^O^Vp 


-1 


and  the  last  term  in  Eq.  (A. 10)  represents  the  effect  of  the  laser  heating 
on  the  change  of  density  in  the  medium.  Since  we  seek  a  wave  equation  for 
the  acoustic  pressure  p,  we  eliminate  u  between  Eq.  (A.l)  and  (A. 2).  It  is 
found  that 


V2p-Ptt+  (— ■ ~  p)y2  pt  =0  •  (A. 11) 

Eliminating  the  p  between  Eq.  (A. 10)  and  Eq.  (A. 11)  yields  a  wave  equation 
for  the  acoustic  pressure  p  and  the  input  entropy  s. 


V2p 


c‘2ptt-stt  po  BJocpl 


(-T) 


c"2pfst  P0 


=  0  .  (A. 12) 


Now  let  us  relate  the  change  in  entropy  s  to  the  heat  q  added  to  the 
medium  by  the  laser  beam.  The  conservation  law  which  governs  that  relation 
is  the  second  law  of  thermodynamics  for  a  reversible  process. 


P  f  ds  =  q  dt 


(A. 13) 


or,  in  linearized  form. 


poVt 


(A. 14) 


Combining  Eq.  (A. 14)  with  Eq.  (A. 12)  to  eliminate  the  entropy  fluctuation  s 
yields 


2  i 

v  p  -  rz  Ptt  * 

c 


X  +2  P 


v2pt =  -^t + 


X  +  2P 


(- 

\Cn 


V^q.  (A. 15) 


This  is  the  Westervelt-Larson  equation  for  a  viscous  (but  not  heat 
conducting)  irrotational  medium  containing  a  heat  source  of  strength  q. 
For  a  lossless  medium,  it  reduces  to  the  classical  form 


V 


c2  P 


tt 


(A. 16) 


47 


*  V  *  T  *  * 


t 


» 


I 


L 


I 


| 


I 


> 


I 


-*  * 

I 


T""  ». 


V  "V*c  » v  -  V - w ■■  Jf ;  J  ^  J  v  *  ”  1  w 


Forming  a  Taylor  series  expansion  of  Eq.  (A. 3)  about  the  equilibrium 
quantities  pQ,  pQ,  and  sQ,  we  have 

p0+  (if-)  (P-po>  *  (ff-)  <*-‘o>  *  (A-4) 

s  p 

By  definition  the  small  signal  sound  speed  is  c  = 

Eq.  (A. 4)  can  be  rewritten  as 

P  =  c_2p  +  s  (— )  .  (A. 5) 

\  3i/p 

The  partial  derivative  in  Eq.  (A. 5)  is  the  result  of  the  change  in  density 
due  to  the  change  in  entropy  induced  by  the  laser  heating.  It  is  therefore 
related  to  the  coefficient  of  thermal  expansion  B  of  the  medium. 

By  definition 


so  that 


or 


(A. 6) 


(A. 7) 


But  it  is  well  known  from  thermodynamics  that  the  specific  heat  at  constant 
pressure  cp  can  be  expressed  as 


c 


P 


(A. 8) 


Combining  Eq.  (A. 8)  with  (A. 7),  and  linearizing  the  result  about  the 
undisturbed  values,  yields 


£p\ 
3s/p 

therefore,  Eq.  (A. 5)  becomes 


-»0 


(A. 9) 
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The  wave  equation  describing  the  pressure  field  of  a  lossless  medium 
containing  a  heat  source  can  be  found  in  Ref.  15.  Starting  from  this 
thermoacoustic  wave  equation,  Westervelt  and  Larson  studied^-  the  special 
case  of  a  laser-induced  sound  field,  discovering  then  the  fascinating 
properties  of  what  is  now  called  a  thermoacoustic  array.  In  this  appendix 
we  rederive  from  the  hydrodynamical  equations  of  motion  the  Westervelt- 
Larson  wave  equation  for  laser-induced  sound,  including  now  the  effect  of 
viscosity. 

Let  P=pQ+p  be  the  total  pressure  in  the  medium,  where  pQ  is  the 
ambient  pressure  and  p  is  the  acoustic  pressure  generated  through  the 
thermal  mechanism.  Similarly  let  p  =  pQ+  p  be  the  total  density  of  the 
medium,  T=TQ+T  be  the  totja  1  temperature,  and  s=sQ+s  be  the  total  entropy. 
The  subscript  o  refers  to  an  undisturbed  quantity. 

Since  the  acoustic  disturbances  generated  by  the  thermal  mechanism  are 
expected  to  be  small,  we  start  directly  from  the  linearized  hydrodynamical 
equations  of  motion:  the  conservation  of  mass  and  the  conservation  of 
linear  momentum. 


mass:  p^  +  p  V  •  u  =  0  ( A . 1 ) 

L  U 

-v  2-* 

momentum:  pQut  +  Vp=(X+2y)Vu  .  (A. 2) 

Here  u  denotes  the  particle  velocity,  (A+2y)  represent  the  viscosity  of  the 
medium,  and  the  subscript  t  denotes  a  time  derivative. 

It  is  implicitly  assumed  in  Eqs.  (A.l)  and  (A. 2)  that  the  medium  is  at 
rest,  steady,  uniform,  and  irrotational ,  and  that  heat  conduction  can  be 
neglected. 

The  thermal  input  from  the  laser  intensity  induces  a  change  in  entropy 
in  the  system,  so  that  an  equation  of  state  describing  the  process  is 

state:  p  =  p(p,s)  (A. 3) 
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(2)  Extension  of  the  theory  to  include  the  effect  of  the  finite 
beamwidth  of  the  laser,  with  a  Gaussian  intensity  distribution  across  a 
horizontal  section  of  the  beam. 

(3)  Generation  of  theoretical  directivity  patterns  in  the  horizontal 
and  vertical  planes. 

Experimental  work  on  moving  thermoacoustic  arrays  is  also  in  progress  and 
will  be  presented  in  a  later  report. 


V.  CONCLUSIONS 


A  time  domain  approach  has  been  developed  for  the  numerical 

computation  of  the  pressure  waveform  radiated  by  a  moving  thermoacoustic 

array.  The  numerical  results  obtained  are  in  good  agreement  with 

g 

analytical  expressions  obtained  from  a  frequency  domain  approach,  valid 
only  in  the  farfield,  and  for  either  very  short  or  very  long  arrays.  The 
major  advantages  of  the  model  described  in  this  report  are  listed  below. 

(1)  The  model  is  valid  in  the  nearfield  of  the  source. 

(2)  It  is  not  restricted  to  the  limiting  cases  of  very  short  or  very 
long  arrays. 

(3)  It  does  not  require  transforms  to  get  time  information. 

(4)  It  is  valid  for  a  source  moving  at  any  velocities  including  the 
transonic  case,  and  it  can  also  be  extended  to  the  case  of  nonuniform 
source  velocities. 

(5)  It  is  not  restricted  in  principle  to  the  case  of  a  source  moving 
rectil inearly,  and  it  seems  particularly  appropriate  in  describing  the 
promising  case  of  a  source  moving  along  any  discontinuous  or  curved 
paths 

Further  work  on  the  theory  presented  here  can  be  divided  into  three 
tasks. 

(1)  Extension  of  the  theoretical  model  to  the  general  case  where  the 
observation  point  is  not  necessarily  in  the  plane  of  motion  of  the  source. 


r"aT  T*  T"  ■*  " 


1  t— -w  ;  v  ‘j  v--  \ 


The  total  pressure  waveform  received  in  the  farfield  can  therefore  be 
expressed  as 


PT(t) «  6  ^t  -  ^  *  I ' (t)  =  I"(t  - 


(50) 


B.  Long  Array 


The  case  of  a  long  and  narrow  thermoacoustic  array  has  also  been 
treated  by  Lyamshev  and  Sedov^  for  farfield  radiation.  They  found  that  the 
pressure  waveform  can  be  written*  in  the  form 

r. 
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where  q  =  ^  I(t)dt  is  the  "area"  of  the  laser  pulse  and  K 

constant  of  proportionality  containing  the  spherical  spreading 

\J 

and  T^cos^/ac.  It  turns  out  that  the  exponential  term  in  Eq.  (51)  is  ot 
negligible  order,  and  keeping  this  term  is  inconsistent  with  some  other 
approximations  assumed  in  their  derivation.  Therefore,  the  pressure 
waveform  can  be  written  simply  as 


P  =  -kls  I 


(52) 


Equation  (52)  describes  an  inverted  laser  pulse.  As  indicated  in 
Fig.  8(b),  the  convolution  approach  described  in  this  report  reduces  almost 
perfectly  to  Eq.  (52)  when  the  source  has  the  shape  of  a  long  and  narrow 
cylinder  and  for  a  receiver  located  in  the  farfield. 


*A  typographical  error  appears  in  Eq.  (13),  Ref.  9,  p.  11.  In  Lyamshev  and 
Sedov's  notation,  the  exponent  of  the  exponential  should  be  divided  by  x  . 
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IV.  COMPARISON  WITH  PREVIOUS  WORK 


The  numerical  results  obtained  by  the  impulse  response  approach 
described  in  the  previous  sections  can  be  compared  with  some  analytical 

Q 

expressions  obtained  by  Lyamshev  and  Sedov  for  the  special  case  of 
farfield  radiation. 

A.  Short  Array 

It  is  shown  in  Ref.  9  that,  in  the  acoustic  farfield,  the  pressure 

★ 

waveform  radiated  by  a  thermoacoustic  array  is  proportional  to  the  second 
time  derivative  of  the  laser  pulse  intensity,  provided  the  effective  length 
(L=  of^)  of  the  array  is  very  small  compared  to  a  typical  acoustic 
wavelength. 

Figure  8(a)  shows  that,  under  these  assumptions,  both  the  numerical 
and  the  analytical  methods  are  in  excellent  agreement.  The  comparison  was 
made  for  a  modulation  frequency  of  5  kHz,  a  laser  pulse  duration  of  1  msec, 
and  for  a  stationary  array.  In  this  case  the  numerical  program  performs 

exactly  a  convolution. 

The  reason  one  should  expect  a  pressure  response  proportional  to  the 
second  derivative  of  the  laser  intensity  can  be  explained  in  the  time 

domain  with  the  impulse  response  approach.  When  the  thermoacoustic  array 
is  stationary,  the  pressure  response  Py(t)  can  be  expressed  as  a 
convolution  between  the  impulse  response  h(t)  and  the  time  derivative  of 
the  laser  intensity.  However,  in  the  geometric  farfield  (L«rQ)  of  a  short 
array  (kL<<  1),  the  impulse  response  tends  to  behave  like  the  time 

derivative  of  a  delta  function  centered  at  time  t=rQ/c.  This  is  clearly 
shown  in  Fig.  4.  The  physical  explanation  for  this  behavior  lies  in  the 
pressure  release  characteristics  of  the  boundary  between  air  and  water. 

*The  "minus"  sign,  on  the  right-hand  side  of  Eq.  (12)  in  Ref.  9,  should  t-e 
read  as  a  plus  sign. 


=  15.6  kHz  at  M=1.5.  This  is  in  reasonable  agreement  with  the  Doppler 
shifted  frequencies  of  13.6  kHz  and  16.0  kHz  in  Figs.  7(b)  and  7(d). 

Finally,  it  is  interesting  to  note  in  Fig.  7(d)  that  time  inversion 
effectively  occurs  at  supersonic  source  velocities.  This  is  due  to  the 
fact  that,  for  a  source  moving  faster  than  the  acoustic  disturbances  it 
generates,  the  first  signal  to  arrive  at  the  receiving  point  will  be  the 
last  one  emitted  by  the  source,  provided  the  source  is  moving  towards  the 
receiver.  This  time  inversion  or  phase  reversal  has  been  known  for  more 
than  a  century.^ 
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